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ABSTRACT 

We estimated the cross-power spectra of a galaxy sample from the Wide-field Infrared Survey 
Explorer (WISE) survey with the 7-year Wilkinson Microwave Anisotropy Probe (WMAP) 
temperature anisotropy maps. A conservatively-selected galaxy sample covers ~13000D°, 
with a median redshift of z=0. 15. Cross-power spectra show correlations between the two data 
sets with no discernible dependence on the WMAP Q, V and W frequency bands. We interpret 
these results in terms of the the Integrated Sachs-Wolfe (ISW) effect: for the |6| > 20° sample 
at I =6-87, we measure the amplitude (normalized to be 1 for vanilla ACDM expectation) of 
the signal to be 3.4± 1.1, i.e., 3.1er detection. We discuss other possibilities, but at face value, 
the detection of the linear ISW effect in a flat universe is caused by large scale decaying 
potentials, a sign of accelerated expansion driven by Dark Energy. 
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1 INTRODUCTION 

Dark Energy (DE) is the principal puzzle of twenty-first century 
physics. When identified with vacuum energy, the quantum field 
theory predic tion for it is either 122 orders of m agnitude too high 
or it is zero dCaldwell & Kamionkowskill2009h . Both predictions 
glaringly contradict astronomical observations suggesting that the 
amount of DE today is of the same order of magnitude as that of 
Dark Matter. The existence of DE rests on measurements of the 
cosmic microwave background fluctuations (CMB), particularly 
by the Wilkinson Microwave Anisotropy Probe satellite (WMAP) 
constraining essentia lly the total energy density in the universe 
dKomatsu et alj|2oTll) . and the low redshift acceleration of th e uni- 
verse detected, e.g., through supernovae dRiess et al. I ll998l) . The 
evidence for DE appears to be overwhelming today, albeit indi- 
rect. Most present and future observations rely on standard candles 
and/or rulers to quantify the geometry and expansion history, or 
growth history of the universe (with the possible exception of grav- 
itational lensing). 



tude smaller than the CMB correlations dCrittenden & Turokl 19961 : 
IPeiris & Spergeill2000h . 

Such correlations were measured in the literature be- 



Luminous Red Galaxies i Scranton et alj 12003; Fosalba et al. 


20031; IPadmanabhan et al.l 20051: iGranett et alj 120081: 


Papai et al. 


201 lb. APM galaxies dFosalba & Gaztaftad 12004 


), infrared 


galaxies l Afshordi et al. 2004b, radio galaxies 


Nolta et al. 



The Integrated Sachs-Wolfe (ISW ISachs & Wolfe] 1 19671 ; 
iRees & Sciama|[T968i) effect promises a unique avenue to directly 
detect the effect of DE on CMB photons: as photons cross potential 
wells, they suffer both blue and redshift cancelling each other in a 
flat, DE-free universe where the potential wells are frozen. How- 
ever, when DE is present, potential wells decay during the pho- 
ton crossing, therefore the redshift fails to fully compensate for 
the blueshift. The net result is a slight "kick" the photon receives. 
This will ultimately result in a tiny correlation of hot spots in the 
CMB with large-scale structure, an effect which is orders of magni- 



jBoughn & Crittenden! l2004blla) . In addition to the ISW effect, 
CMB photons become correlated with the LSS at high I (small 
scales) due to the Sunya ev-Zeldovich effect , and lensing. The latter 
was recently detected bv lSmithetai] ( l2007b . 

This analysis uses the fir st data release of th e Wide-field In- 
frared Survey Explorer (WISE: lwright et ap2010b . We use a sam- 
ple of galaxies over a survey area of ~13000 D° with a median 
redshift of z=0. 15. The combination of area and volume at mod- 
erate redshift has the potential to detect (or reject) the ISW effect 
more clearly than previous studies. 



2 ANALYSIS AND RESULTS 
2.1 Data 



For th e CMB, we use the WMAP 7-year data set dKomatsu et alj 
|20T 2). We use the Q, V and W foreground removed maps. We re- 
bin the m aps in Healpix form at to n s idc = 128 for a pixel size of 
27'/pixel jGorski et alj|2005b . We apply the Temperature Analysis 
Mask provided by the WMAP team to exclude regions of Galactic 
emission and known point sources leaving about 78% of sky. 
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Figure 2. Galaxy number density (in 27'pixels) as a function of Galactic latitude b and ecliptic coordinates /3 and A. The red line connects median values, 
with RMS in the red vertical lines. 



Galaxy Counts 




Figure 1. Galaxy counts from the WISE galaxy sample. Regions with zero 
galaxy counts are masked out by the object masks used for the WISE galaxy 
and WMAP sample. These are regions with |6| < 10 deg, missing data, and 
star clusters. 
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Figure 3. The normalized redshift distributions of the WISE galaxy sample 
(red solid line) and the GAMA spectroscopic sample (black dotted line). 



The WISE satellite surveys the entire sky at wavelengths of 
3.4, 4.6, 12, and 22 /im (Wl through W4), with the point source 
sensitivities of 0.08, 0.11, 1, and 6 mjy, respectively. We use the 
first data release issued on April 14, 2011, containing about half 
the sky. 

We need to remove stars from the WISE object catalogue. We 
find that most stars (|6| < 10°) have a 3.4 — 4.6/im color less than 
0.2. Therefore, we apply a color cut of 3.4 — 4.6jj,m >0.2 to select 
galaxies. We further apply a color selection of 4.6 — 12/im >2.9, 
which removes a greater number of stars at the expense of a small 
number of galaxies. Due to the sun-synchronous orbit of the satel- 
lite, WISE data are deeper in pole regions. To select a uniform 
galaxy sample, we only use galaxies with Wls.e^m < 15.2 mag. 

We construct Healpix maps of the galaxy counts matching 
our CMB maps, shown on Fig. [T] We apply the WMAP mask 
and additionally mask out regions with missing data and unusually 



high object-counts such as star clusters and the Magellanic clouds. 
These can be seen as thin stripes in Fig.Q] totaling 613 0° of area. 
We also mask cells with more than 10 flagged sources (as persis- 
tent, halo, ghost, or spike). This helps to exclude cells contami- 
nated by bright nearby galaxies and stars, as visual inspection of 
pixels in the tail confirmed. We do not use regions with |6| < 10° 
due to strong stellar contamination. Our galaxy catalog contains 1.2 
million galaxies over an area of 13622 D° and the median galaxy 
density is 86 deg -2 . 

We investigate gradients in the galaxy density in Fig. [2] as a 
function of Galactic latitude 6, and ecliptic coordinates /? and A. 
There appears to be a significant decrease in the galaxy number 
density at 10 < \b\ < 15, reaching 30% at |6| = 10°. At this 
low galactic latitude, Milky Way stars can mask background galax- 
ies and co nfuse photometry . This effect has been studied in opti- 
cal bands (Ross etal .1120 111) , but it is more severe in WISE due to 
the broader point spread function. In the middle panel, the galactic 
centre is at A ~270°, and the same effect can be seen. To test how 
this might affect our results, we analyse maps with three different 
galactic cuts (|6| >20°, 15°and 10°) in Section|2~3l The WISE data 
become deeper toward the pole regions, but there is no significant 
change in galaxy density as a function of /3 (right panel). 

To estimate the expected ISW signal, we need to know the 
redshift distribution of our galaxy sample, and the galaxy bias. We 
have cross-matched t he WISE galaxy sa mple with GAMA (Galaxy 
and Mass Assembly; iDriver et alj|201 lh sample to estimate a pre- 
liminary redshift distribution. GAMA is a spectroscopic sample of 
205,000 galaxies observed with AAf2 4m telescope. The survey 
depth is as deep as tab < 20.5 mag, and thus, deep enough to es- 
timate the redshift distribution of the WISE galaxy sample. In the 
overlapping regions, 83% of the WISE galaxy sample had spectro- 
scopic counterpart in GAMA. Fig.[3]shows the redshift distribution 
of the WISE galaxy sample with GAMA counterparts. The dot- 
ted black histogram shows the redshift distribution of the GAMA 
sample, which has a much longer tail toward higher redshift. The 
median redshift of the WISE-GAMA sample is z=0.148. 

2.2 Galaxy bias 

In order to determine galaxy bias, we measure the galaxy-galaxy 
angular power spectrum. The form of the angular power spec- 
trum can be theoretically computed given cosmological parame- 
ters. We use the WMAP 7-year ACDM cosmological parameters 
(Ho = 70.4 km/s/Mpc, fi m = 0.272). However, in linear theory, 
bias cannot be determined independently of erg, the power spec- 
trum is normalized as Ci oc (bag) 2 , where b relates galaxy and 
matter overdensities as 5 g = bS m - We therefore use a fixed value 
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Figure 4. Auto correlation of the WISE galaxy sample (|6| > 20°). The 
red solid curve is a fit of a theoretical matter power spectrum from CAMB 
with the latest cosmological parameters for ACDM. The gray solid line is a 
Poisson shot noise. 



of crg=0.80, and fit the galaxy-gala xy angular correla t ion only for 
bias. We use CAMB with Halofit jLewis et alJbood : ISmith et all 
120031) to generate non-linear matter powerspectra for our cosmo- 
logical model, and obtain three-dimensional spectra at the median 
redshift of z=0.148. The angular power spectrum is a projection of 
the three-dimensional power spectrum with a kernel that depends 
on the redshift distribution: 



' U TV 



<t>(r)ji(kr)dr 



k 2 P(k)dk, 



(1) 



with cj>(r) oc dN ^ r > jjZL for comoving coordinate r normalized such 
that J cf>(r)r 2 dr — 1 and ji is the spherical Bessel function. We 
make use of the small-angle limit with the L imber equation to eva l- 
uate the power spectrum at / > 20, see e.g., lAfshordietaT]j2004l) . 

We compute the galaxy-galaxy angular power spectrum of the 
WISE galaxy sa mple using SpICE(Sp atially Inhomogeneous Cor- 
relator Estimaton lSzapudi et al.l2001al lbT). a fast quadratic estimator 
which inverts the coupling matrix of the Cj 's in pixel space where 
it is diagonal. The resulting power spectrum of the WISE galaxy 
sample is shown in Fig. [4] The Poisson shot noise is given by 1/N, 
where N is the mean number of galaxies per steradian. In our sam- 
ple, 1/N =3.4xl0 -6 , which is much smaller than the correlation 
at I = 2 — 100. We subtract the Poisson noise (gray line), then fit 
the amplitude of the model template over the range 1 = 2 — 100. 
We find a bias parameter of b 9 = bag/O.S = 1.06 ± 0.05. 



2.3 WISE- WMAP Cross Power Spectrum 

We cross-correlate the WISE galaxy sample with the foreground- 
removed Q, V, and W maps from WMAP. Because of the declining 
galaxy density at low galactic lattitude b (Fig|2j, we show results in 
three different samples, |6| > 10, |6| > 15 and |6| > 20. But first, 
we show results with the |6| > 20 sample. The resulting angular 
power spectra are shown in Fig. [5] The spectrum (black dots) band 
powers are binned 0.15 dex logarithmic bins. The boundaries are 
Z=6,8,l 1,16,22,31,44,61 and 87 such that the first band includes 
I =6,7, etc. We do not use data at I ^87 since these small scales 
should not be sensitive to the ISW effect. 

In order to estimate uncertainty in the WISE- WMAP cross 
power spectrum, we have run 1000 simulations. Using synf ast 
and the latest ACDM cosmological parameters, we create 1000 ran- 
dom CMB maps. Cross-correlating these random CMB maps with 



the WISE galaxy catalog, we estimated the covariance matrix. It 
is dominated by diagonal terms, i.e. neighboring band powers are 
uncorrelated to a good approximation. 

The error bars plotted in Fig.[5]are given by the diagonal ele- 
ments of the covariance matrix. The large-scale modes at / < 6 are 
not well constrained on the cut sky. The size of a contiguous sur- 
vey region in the WISE map (Fig.[T} is ~60°, or 30° in radius. This 
corresponds to a limiting scale of I ~ 180° /9 ~6. Therefore, we 
do not have strong constrains at I < 6. For simplicity we exclude 
the large-scale modes at I < 6 in the analysis. 

The C'e's are consistent among different WMAP frequencies, 
while most foregrounds from Milky Way emission should exhibit 
some color dependence. 

Theoretical expectation for the ISW power spectrum from the 
the latest cosmological parameters for ACDM are plotted with the 
red solid line. We adopted bias the parameter of b 3 = 1.06 es- 
timated from auto-correlation in Section 12.21 The galaxy density- 
CMB cross power spectrum is given by (Cooray 2002), 



Cf = Tcm, 
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k dkP k x 



r' 2 (f>(r')ji(kr')dr' (2) 



where Di(z) is the linear growth factor. 

The measured power spectrum is higher than the ACDM pre- 
diction. To quantify this we scale the model with a free amplitude 
parameter. In Fig. [5] the best x^-fA was in the orange line, which is 
2.2 a larger than the original theoretical expectation. 



2.4 Significance tests 

We now test the following three hypotheses for consistency with 
the data: (a) theoretical best-fit ISW (orange line in Fig. [5j, (b) 
theoretical ISW from vanilla ACDM (red line) (c) Null hypothesis: 

no power, i.e. n o ISW effect. 

We follow iFrancis & Peacockl d2010h and compute \ 2 values 
using the covariance matrix for the three hypotheses. In Fig. [5] we 
see a high power at 1=9 and 12, that drops at I — 7. The signifi- 
cances of correlation at I — 9 and 12 are 1.3 and 1.7cr, respec- 
tively. While the significance of an individual band power is af- 
fected by the binning, the final x 2 values and likelihoods are robust. 
To check for this, tested alternative band powers with boundaries at 
1=6,9, 12, 17,24,48,67 and 95 (purple squares in FigO; none of our 
numerical results or conculsions changed. 



The \ 2 statistic is given by x 2 



(x - t i ) T C- 1 (x - ti) 
for the observed data vector x and a given theoretical expectation 
tj for model i. For Z=6-87 with 8 data points, the \ 2 values are 
11.6, 16.1 and 20.7 for (a), (b), and (c), respectively. The \ 2 values 
decrease significantly for the ISW models (a) and (b) compared 
with the null hypothesis. Computed Ax 2 (Ratio of evidence) are 
9.2 for Z=6-87. Generally, Ax 2 ^3 can be considered as strong 
evidence. Our best-fit measurement has Ax 2 >3 in all cases. 

In terms of the likelihood ratio, e 1//2Ax , we obtain 97 and 
9.9 as shown in Tabl^T] The likelihood ratio is considered signifi- 
cant if it is greater than 5. According to these ratios of evidence, the 
ISW hypotheses (a) and (b) are preferred over the null hypothesis. 

For the model (a), these numbers can be compared with the 
significance of the amplitude fit itself. We computed 1-er errors of 
the fit of the theoretical ISW model to the data, by estimating the 
68 percentile of the e _1,/2x distribution. We also checked that the 
measured errors are consistent with the width of the best-fit Gaus- 
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sian. Results in Table[2]show that the best-fit model is 3.1a from 
the null hypothesis. 

2.5 Systematic tests 

We have shown results with |6| > 20° sample. In this section, we 
explore if changing the galactic cut affects our results. 

With a |6| > 15° cut, we have ~1700 D° more area than the 
|6| > 20° cut. The resulting amplitude of the best-fit ISW model is 
3.2±1.1, i.e., 2.9 a, in agreement with the result from |6| > 20°. 
As shown in Table Q] the A% 2 and likelihood ratio are decreased, 
but are not significantly changed. At \b\ ~ 15° the galaxy density 
starts to decline towards the galactic plane, and thus, the correlation 
signals have been reduced some, but overall, we obtained consistent 
results with the |6| > 20° analysis. 

Next we used |6| > 10° sample. This sample adds another 
~ 1600 n°. As shown in Tables [T]and[2] the signal becomes lower 
and the error is larger giving a worse significance of 2.0<r. This is 
most-likely due to the systematic decrease in galaxy density at low 
b, where artefacts of bright stars and high stellar density decrease 
the galaxy density. Large-scale gradients in the galaxy distribution 
can increase the measurement errors because they are amplified by 
cosmic variance of the CMB on large angular scales. 

To test this interpretation, we artificially flattened the galaxy 
density as a function of b, by dividing the galaxy density by the 
mean galaxy density at each 6. This process can be considered as 
a high pass filtering that removes large-scale gradients. As shown 
in Tables Q] and [2] the significance is increased to 3.2 a, becom- 
ing consistent with values we measured with |6| > 20° and 15° 
samples. 

To the opposite direction, we tried |6| > 25°, again obtain- 
ing a consistent amplitude of 3.5±1.2 (Tables[T]and[2]l. These tests 
show that our results are robust, obtaining ~3 a significance re- 
gardless the choice in the galactic cut b. 

We briefly discuss stellar contamination next. Stars are not 
clustered, but star counts might be correlated with the CMB due 
to galactic contamination. According to Fig. [2] the galaxy density 
is constant at |b| >20, although there is a slight negative gradi- 
ent at |fo| <20. We argue that is likely to be due to artefacts and 
confusion with stars. The robustness of the results with respect to 
galactic cuts, however, implies that the effect on the correlations is 
negligible for the foreground subtracted CMB maps we use. We do 
not see any significant colour dependence either between WMAP 
frequencies (V, W, and Q) in Fig. [5] These points suggest that the 
stellar contamination is small, or at least uniform across our sur- 
vey area. In the worst case the contamination would be 17%, if all 
unidentified GAMA sources were stars, but it is likely to be much 
less. 
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Figure 5. Power spectra between the WISE galaxy sample (|fe| > 20°) 
and WMAP (V,Q,W) maps. The data in original resolution (black dots) are 
binned in 0. 15 dex logarithmic bins and shown with error bars. The red solid 
line is theoretical expectation with the bias of 1.06, while the orange line is 
the best fit to the observed data. The fit used 8 (7=6-87) data points. The 
purple line is from WMAP W data in different binning. 



Table 1. Significances from x teste at I = 6 — 87. 



Galactic cut 




x 2 


d.o.f. 


Ax 2 


Likelihood ratio 


\b\ > 25° 


Best-fit ISW 


10.2 


7 


8.7 


77 




ISW with ACDM 


14.7 


8 


4.2 


8.2 




Null hypothesis 


18.9 


8 







\b\ > 20° 


Best-fit ISW 


11.6 


7 


9.2 


97 




ISW with ACDM 


16.1 


8 


4.6 


9.9 




Null hypothesis 


20.7 


8 







|6] > 15° 


Best-fit ISW 


11.0 


7 


8.1 


59 




ISW with ACDM 


14.9 


8 


4.3 


8.6 




Null hypothesis 


19.2 


8 







|fe| > 10° 


Best-fit ISW 


9.1 


7 


3.8 


6.6 




ISW with ACDM 


11.1 


8 


2.6 


3.2 




Null hypothesis 


13.7 


8 







\b\ > 10° 


Best-fit ISW 


14.3 


7 


10.2 


166 


with flattening 


ISW with ACDM 


19.1 


8 


5.4 


14.8 




Null hypothesis 


24.5 


8 








3 DISCUSSION AND CONCLUSIONS 

The significant correlations between WMAP and WISE on large 
scales may have contributions from many sources including the 
Milky Way, extragalactic point sources, even zodiacal light. We 
have attempted to minimize contamination by constructing a clean 
galaxy sample from WISE. We find no systematic trends in galaxy 
number density with ecliptic latitude and longitude (Fig. |2j. We 
use the WMAP foreground reduced CMB maps to reduce the sen- 
sitivity to Milky Way emission in WMAP bands. The primary 
sources (dust, synchrotron and free-free emission) have character- 
istic spectral shapes across the WMAP frequency bands. The ex- 



Table 2. Amplitude and errors of the best-fit ISW model at £=6-87. 



Sample 


Area (tJ°) 


Amplitude and error 


a 


|6| > 25° 


8495 


3.5±1.2 


3.1 


6 > 20° 


10337 


3.4±1.1 


3.1 


6 > 15° 


12032 


3.2±1.1 


2.9 


|6 > 10° 


13622 


2.6±1.3 


2.0 


|6| > 10° with flattening 


13622 


3.2±1.0 


3.2 
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cellent agreement between the Q, V and W spectra (Fig. [5} leads us 
to believe that Galactic foregrounds are unimportant. 

Extragalactic point sources below the WMAP detection limit 
could also contaminate the signal. Since WISE is a near-IR sur- 
vey, galaxies detected in WISE may also emit strongly at WMAP 
frequencies. Such sources would produce a characteristic signal in 
the power spectrum that could become important on small scales. 
The zeroth order estimate of this effect is a Poisson term with a 
fiat C; spectrum; no evidence for this is seen at / < 87. Simi- 
larly, the Sunyaev-Zeldovich effect is expected to appear at higher 
I and has negligible effect at I < 87 dFosalba & Gaztanadl2004l ; 
lAfshordi et al]2004) . 

We leave detailed discussion of the cosmological significance 
to a future work, a few comments will suffice here. Our best fit ISW 
model is 2.2 a higher than ACDM expectation. This i s consis- 
tent with many previous cross-correlation measure ments jHo et al.l 
l2008l : lGiannantonio et al.l2008MGranett et alj2009h that had higher 
amplitude than vanilla ACDM prediction. Possible reasons include 
uncertainties in cosmological parameters or our galaxy bias model. 
Additionally, the measurement is subject to significant cosmic vari- 
ance originating from both the CMB and LSS fields jPapai et all 
1). At present there is at most a mild tension between vanilla 
ACDM and measurements that is not fully understood. 

In the ACDM model, the amplitude of the ISW signal is well 
constrained by our knowledge of erg, the Hubble constant and the 
matter density. However, we do face uncertainty in the redshift dis- 
tribution and its effect on our estimate of the linear bias and the 
amplitude of the ISW spectrum. We obtained the dN/ dz by match- 
ing the WISE galaxy sample to the GAMA spectroscopic sample. 
Although the GAMA sample is deep, it is an optically selected sam- 
ple while WISE is sensitive to the near IR. Thus, GAMA may be 
missing a population of WISE galaxies such as heavily obscured 
one. We checked the fraction of our WISE galaxy sample detected 
by GAMA was 83%. While 83% is a large fraction, the redshift 
distribution of the rest 17% is unknown. 

However, even if we assume this 17% of galaxies are all at 
higher redshift than 0.148, the median only changes to z=0. 155. 
More generally, to test whether erroneous redshift distribution can 
affect bias estimate, we have artificially shifted the median red- 
shift from z=0.09 to z=0.22, a value clearly beyond reasonable. The 
change in bias is within 40%. Therefore, the errors in the redshift 
estimation alone can only account for a fraction of the shift between 
the original and best-fit theoretical ISW estimations. Furthermore, 
the amplitude of the ISW spectrum is relatively insensitive to the 
median redshift over this range. 

The most relevant measurements in the literature to com- 
pare with are from 2MASS for which we ak and null detec- 
tions of the ISW signal have been claimed jRassat et alj |2007| ; 
iFrancis & Peacockj201fj|) . However, the median redshift of 2MASS 
is z=0.07 and in terms of comoving volum e, the WISE sample is 
4.6 times larger. IFrancis & Peacock! d2010h selected a higher red- 
shift sample (0.2 < 2 < 0.3) using photometric redshifts but it is 
also much more sparse. The WISE sample may be situated in an 
ISW 'sweet spot': it is a large volume and samples a low redshift 
when dark energy has come to dominate the cosmic expansion. In- 
deed, the WIS E survey seem s in better agreement with higher red- 
shift datasets dHo et al ] |2008l : iGiannantonio et alj|2008h . However, 
the ACDM expectations are at the detection limit for current sur- 
veys. It is possible that we are only measuring the ISW signal from 
"lucky" statistical fluctuations that push the signal above the typi- 
cal 2-cr detection threshold leading to biased conclusions if taken 
at face value. Future surveys that expand on the sky coverage and 



redshift range in the next decade will provide a clearer picture on 
the integrated Sachs-Wolfe effect and the role of dark energy. 

We thank the anonymous referee for many insightful com- 
ments. We acknowledge financial support from the NASA grant 
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Office for the Research and Technology (NKTH). 
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